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1 „  Introduction 


Recently,  a  continuum  theory  of  finitely  deformable,  heat  conducting 

composite  materials  was  developed  by  modeling  N  identifiable  constituents 

•* 

as  interpenetrating  solid  continua-.  In  deriving  the  general  system  of 

nonlinear  equations  governing  the  behavior,  the  motion  of  a  point  of  the 

combined  continuum  was  permitted  to  be  finite  while  the  relative  motion  of 

the  individual  constituents  was  constrained  to  be  infinitesimal  in  order 

that  the  solid  composite  not  rupture.  The  restriction  imposed  in  Ref.l  which 

demands  that  the  relative  motion  of  the  constituents  be  infinitesimal  is 

2  _s 

one  of  the  features  that  distinguishes  this  theory  from  other  work  on 

composites.  Another  important  feature  distinguishing  the  description  in 

2-8 

Ref.l  from  that  of  tha  other  work  is  that  r.o  energy  of  interaction  be¬ 
tween  the  constituents  is  omitted  in  Ref.l  as  it  is  in  the  other  work. 

In  this  payer  we  specialize  the  theory  developed  in  Ref.l  and  consider 
the  one-dimensional  motion  of  a  two  constituent  composite  material  whose 
identifiable  constituents  are  elastic,  as  is  the  interaction  between  the 
constituents  with  the  exception  of  the  volumetric  part.  We  examine  the 
behavior  of  one-dimensional  acceleration  waves  in  such  media  on  the  assump¬ 
tion  that  thermodynamic  influences  may  be  ignored.  In  particular  we  seek 
to  determine  how  the  behavior  of  such  waves  is  influenced  by  (i)  the  mechan¬ 
ical  properties  of  the  mean  (center  of  mass)  behavior  of  the  combined  com¬ 
posite,  (ii)  the  mechanical  properties  associated  with  the  relative  motion 
of  the  individual  constituents,  (iii)  the  coupling  between  these  two  motions, 
(iv)  the  relative  mass  densities  of  the  individual  components  of  the  composite, 
and,  finally,  (v)  the  dynamical  conditions  prevailing  ahead  of  the  wavefront; 
and  in  certain  interesting  simplified  special  cases  by  (vi)  the  mechanical 
properties  of  the  individual  constituents  of  the  composite  and  (vii)  the 
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coupling  effects  arising  from  the  volumetric  interaction  between  the  elements 
of  the  composite. 

Section  2  of  this  paper  is  devoted  to  a  brief  review  of  the  equations 
which  govern  the  one-dimonsional  motion  of  two  constituent  composites. 

After  recording  the  global  forms  of  the  equations  which  govern  the  balance 
of  linear  momentum,  we  state  the  constitutive  equations  which  govern  the 
one-dimensional  motions  of  elastic  composites  which  are  made  up  of  two 
identifiable  clastic  media.  Sections  3  and  4  are  devoted  to  the  study  of 
the  propagation  of  acceleration  waves.  In  Section  3  we  show  that  the 
balance  laws  and  constitutive  equations  set  forth  in  Section  2  imply  the 
existence,  in  general,  of  two  distinct  types  of  acceleration  waves,  the 
fast  one  of  which  is  associated  with  the  mean  elasticity  of  the  combined 
composite  and  the  slow  one  with  the  elasticity  associated  with  the  relative 
motion  of  the  individual  constituents.  When  the  effects  of  coupling 
between  the  center  of  mass  motion  of  the  combined  composite  and  the  rela¬ 
tive  motion  of  the  constituents  is  small,  one  wave  propagates  with  a  velocity 
which  is  close  to  that  of  the  ordinary  elastic  wave  speed  of  the  combined 
composite  while  the  speed  of  propagation  of  the  second  wave  is  close  to 
that  of  the  wave  of  the  relative  motion  of  the  constituents.  The  behavior 
of  waves  in  a  number  of  highly  restrictive  special  types  of  composite  is 
examined  in  Section  4.  It  is  shown,  in  particular,  that  when  one  of  the 
components  is  a  chopped  fiber,  only  one  acceleration  wave  may  exist  in  the 
composite  and  when  the  interaction  between  the  constituents  is  purely 
volumetric,  the  two  acceleration  waves  propagate  with  the  respective  speeds 
of  those  in  the  individual  constituents. 

The  manner  in  which  the  amplitudes  of  acceleration  waves  vary  as  they 
traverse  the  composite  is  examined  in  Sections  5  and  6.  A  standard  analysis 
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is  employed  in  Section  5  to  show  that  the  amplitude  of  an  acceleration  wave 
satisfies  an  equation  of  Bernoulli  type.  The  various  typos  of  behavior 
possible  in  a  number  of  situations,  including  the  possibility  of  shock 
formation,  are  discussed  in  Section  6.  In  tho  general  case  where  two  waves 
may  exist  in  the  material  the  behavior  of  both  the  "fast"  and  "slow"  waves 
is  discussed.  It  is  noted  that  while  the  medium  ahead  of  the  fast  wave  may 
be  in  a  steady  state  before  the  arrival  of  the  wave,  this  condition  will 
be  unlikely  to  prevail  ahead  of  the.  "slow"  wave  because  of  the  motion 
induced  ahead  of  this  wave  by  the  passage  of  the  precursor.  The  propaga¬ 
tion  of  a  "fast"  wave  in  a  composite  which  is  initially  at  rest  in  an 
arbitrary  permissinle  state  of  deformation  is  examined  in  detail  and  it  is 
shown  that  the  behavior  of  the  amplitude  of  such  a  wave  is  the  same  as  that 
of  an  acceleration  wave  propagating  in  a  single  phase  elastic  material  which 
is  in  a  state  of  nonhomogeneous  deformation  ahead  of  the  wave,  A  similar 
situation  prevails  when  the  center  of  mass  deformation  is  homogeneous  but 
the  deformation  fields  of  the  two  continua  which  make  up  the  composite 
are  not.  If  the  material  ahead  of  the  wave  is  in  its  natural  stress-free 
state  then  it  is  found  that,  as  far  as  " fast"  waves  are  concerned,  the 
material  behaves  in  the  same  way  as  would  a  single  phase  thermoelastic 
medium. 

The  behavior  of  the  amplitudes  of  acceleration  waves  in  the  highly 
special  cases  treated  in  Section  4, as  well  as  the  higher  order  discontin¬ 
uities  induced  by  some  acceleration  waves, is  also  examined  in  Section  6, 
Thus,  when  the  composite  is  such  that  the  interaction  between  the  consti¬ 
tuents  is  purely  volumetric  arj  depends  only  on  the  relative  displacement 
of  the  constituents,  the  composite  behaves,  as  far  as  acceleration  waves 
arc  concerned,  as  an  elastic  material  composed  solely  of  one  of  the 
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component  continua,  In  this  particular  situation,  the  acceleration  of 
only  one  of  the  components  suffers  a  discontinuity  at  either  of  the 
waves  which  may  exist  and  the  acceleration  of  the  points  of  the  second 
continuum  are  continuous  along  with  the  first  order  time  derivative,  but 
the  second  order  time  derivative  of  the  acceleration  of  the  second  con¬ 
tinuum  suffers  a  jump  discontinuity  at  the  wavefront.  Finally,  we  examine 
the  behavior  of  an  acceleration  wave  propagating  in  a  composite  one  of 
whose  component  continua  is  made  up  of  chopped  fibers.  Only  one  accel¬ 
eration  wave  may  exist  in  such  a  material  and  the  acceleration  of  the 
chopped  fiber  continuum  is  continuous  at  the  wavefront.  The  order  of 
the  discontinuity  in  the  motion  of  the  chopped  fiber  continuum  depends 
on  the  nature  of  the  composite  and  the  conditions  prevailing  ahead  of 
the  wave.  In  general,  the  first  derivative  of  the  acceleration  of  the 
chopped  fiber  continuum  suffers  a  jump  discontinuity  across  the  wave- 
front.  On  the  other  hand,  if  the  composite  is  centrosymmetric  and  is  in 
a  state  of  equilibrium  ahead  of  the  wave  then  the  first  derivative  of 
the  acceleration  of  the  chopped  fiber  continuum  is  continuous  everywhere 
and  for  all  time  while  the  second  derivative  suffers  a  jump  discontinuity 
at  the  acceleration  wave. 

2 .  Basic  Equations  for  One-Dimensional  Motions 

We  are  interested  here  in  studying  the  motion  in  one  dimension  of 
a  composite  consisting  of  two  interpenetrating  solid  continua.  Initially, 
the  two  continua  occupy  the  same  region  of  space  and  hence  the  location 
of  the  identifiable  components  of  the  composite  may  be  specified  by  a 
single  reference  coordinate  X.  It  should  be  noted  that  X  specifies  the 
position  of  a  point  of  each  of  the  interpenetrating  continua  at  some 
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fixed  time  t =  0,  say.  The  subsequent  motion  of  the  composite  is  described 
by  specifying  two  functions 

y(1)  =y(1>  <X,t),  y(2)  =y(2'  <X,t)  ,  (2.1) 

which  give  the  positions  at  time  t  of  the  points  of  the  two  interpenetrating 
continua  which  were  simultaneously  located  at  the  point  X  at  time  t=0. 

V)a  denote  the  mass  density  of  the  ith  constituent  in  the  reference  con¬ 
figuration  by  0^^  and  in  the  current  configuration  by  p^.  The  center 
of  mass,  at  time  t,  of  the  particles  of  the  continua  which  simultaneously 
occupied  the  point  X  at  time  t=0,is  given  by 


y  =  y(X,t)  - 


(1)va)(X.t)+p-2)y(2)(X.t) 


(1)  .  (2) 
P  +  P 


Clearly,  Eqs.  (2.1)  may  be  written  in  the  form 


=  y  +  w(1*  (X,  t) ,  y  <2)  =  y  +  w(2)  (X,  t)  ,  (2,3) 

where  w'"^  (X,t)  is  the  displacement  of  the  point  X  of  the  ith  continuum 
relative  to  the  center  of  mass  of  the  points  originally  at  X  at  time  t=0. 
Rs  in  Ref.l,  we  place  no  restriction  on  the  magnitude  of  y,  but  the  rela¬ 
tive  displacements  w^,  »/2^  are  taken  to  be  infinitesimal.  The  deforma¬ 
tion  gradients  at  the  point  X  are 

F  =  F(X,t)  =aj£y(X,t),  F(l)  =F(1)  (X,t)  =  axy(l)  (X,t)  =F  +  F(;L)  (2.4) 


where 


F(l)  =3xw(l)  (X,t)  . 


In  (2.4)  F  is  the  deformation  gradient  of  the  center  of  mass,  F  is  the 
deformation  gradient  of  the  point  of  the  ith  constituent  which  was  located 


at  the  point  X  at  t=0  and  Fv  is  the  relative  deformation  gradient  of 
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this  point.  Since  w^(X,t)  is  an  infinitesimal  displacement  field, 

| p ( 1 > |  «  |  F |  and,  since  mass  is  conserved  separately  for  each  constituent, 
we  have 

p^l}  =p(l)  (F  +  F(l>)  ~pU)F,  (2.6) 


so  that  we  may  write 

(1)  ^  (2^  ~  ,,  *7 \ 

P  =  P  +  P  ,  pF=  po  ,  (2.7) 

(1)  A  (2) 
p  =  p  +  p  , 

where  p  is  the  total  reference  mass  density  of  the  composite, 
o 

Since  y(X,t)  is  the  position  at  time  t  of  the  center  of  mass  of  the 
points  of  the  constituent  continua  which  were  at.  X  at  t=0,  it  follows 
from  (2.2),  (2.3),  (2.6)  and  (2.7)  that 


rw(1)+w(2)=0 


(?-.8) 


/ i)  (2 )  9 

where  r  =  p'  /o  is  assumed  to  be  constant  .  At  this  point  it  should 
yo  '  o 

be  noted  that  p^^  and  do  not  represent  the  actual  mass  densities 

of  each  of  the  constituents  in  the  composite,  but  only  represent  those 

quantities  in  each  of  the  interpenetrating  continua,  which  occupy  the 

same  region  ot  space  and,  respectively,  represent  each  constituent  in  the 

model.  Suppose  that  at  time  t=0,  the  ith  constituent  occupies  a  fraction 

X.  of  the  volume  of  the  constituent  so  that  p^  where  is 

1  rO  1  o  7  *o 

the  mass  density  which  a  body  composed  solely  of  the  ith  constituent 
would  have.  It  follows  from  (2,7)  that 


-d)  -(2)  ,  .  .  . 

po  =  Xl°o  X2°o  >  Xl  +  X2=1’ 


(2.9) 


while  r  =  RX^/X2,  where  R=  represents  the  constant  ratio  of  the 

actual  mass  densities  of  the  constituents. 
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The  one-dimensional  version  of  the  integral  forms  of  the  equations 
of  balance  of  the  composite  follows  from  Eqs.  (6.1)  and  (6.2)  of  Ref.l  in 
the  form 


dt 


d  r 

at  J 


O' 


(1) 


J  porw'^'d>:  =  £(X  t)  +  J  3dX  , 

x  x 

a  a 


(2.10) 


(2.11) 


where  X^,  are  two  arbitrary  points  in  the  reference  configuration  of 
the  composite,  and 


(1)  .  (2) 


K  =  T  +  T 


ft  (1)  .P  (2) 

Jr  =  t  -  r  r 


3=  <1+  r)  V2 


(2.12) 

(2.13) 

(2.14) 


In  (2.12)  -  (2.14)  K,3  represent  the  total  stress  and  the  relative  stress 

for  the  combined  continuum,  respectively,  and  are  the  stresses 

L  12 

for  each  of  the  interpenetrating  continua,  while  F  is  the  force 
exerted  by  continuum  2  on  continuum  1.  In  Eqs.  (2.10),,  (2,11)  and  in  what 
follows  a  superposed  dot  denotes  material  differentiation:  G  =  cMl  (X,  t)/?it. 

In  addition  to  the  foregoing  we  have  the  relevant  constitutive 
equations’^,  which  we  take  in  the  form 


k=k(f,f(1,j  w(1)),  ■0  =  £(f,f(1),  w(1)) 


3  =  3  (F,  F(1) ,  W(1))+3(F,F'^,  wVX';  v/U )  , 


(1)  (1)  .  (1) 


(2.15) 


.  «  i  .  a 

and  we  assume  that  the  functions  K  (.,.,.),  Jj 3  (.,.,.)  and  3 (.,.,.;.) 
(2) 

are  C  functions  of  their  arguments.  For  future  reference,  we  note  that 
it  has  beer,  shown  in  Ref.l  that  K,  Jt  and  3  are  related  to  the  stored 
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-V  M)  (1) 

energy  density  E  =  E( F,  F  ,  w  )by  the  formulae 


K=po^(P,F(1),w(1)),  ^=Podp(l)S(F,P(I),W(1))  , 

J»-pS  (l)E(F,F(1),w(1))  . 
o  w  ’ 


(2.16) 


For  a  positive  rate  of  entropy  production  Vf ;w  )  must  be  an  odd 


_  4  .  „  . (1)  . 
function  or  w  ,  i.e., 


.(1),  ft.,  .(1). 

3  (.,.,.  ?w  )= -3 w  )  , 


(2.m 


and  the  condition  (2.17)  implies  that  3(-)  has  the  representation 


£.  *  ,  .  (1)  .  (1) 

3  =  g(., .  ;w  )w  , 


(2.18) 


where  g  is  an  even  function  of  w(1^  which  must  be  strictly  negative1, 

Und^r  certain  circumstances  it  turns  out.  to  be  convenient  for  interpretive 
purposes  to  take  the  stored  energy  E  in  the  equivalent  form  E=E(F^, 

\  wa)),  chen  it  follows  from  Eqs.  (2.4),  (2.8),  (2.12),  (2.13)  and 

(2.16)  that 


Tl=Po'di(U)E,  T2-po^r(2)E 


(2.19) 


it  is  clear  that  the  stored  energy  density  may  be  written  in  the  form 


3oE-x1p'1)E<1)(i(1))  +  \2pf)E(2)(F(2): 

+  X1X2poE(12)  (F(1),F(2)  ;wa))  , 


(2.20) 


where  E^1' (F^)  is  the  energy  density  which  the  body  would  have  at  the 
point  X  if  it  were  composed  solely  of  the  ith  constituent.  Wie  third 
term  on  the  right-hand  side  of  (2.20)  may  be  called  the  interaction  energy 
density  of  the  constituent  continua  and  it  is  the  presence  of  this  term 
which  causes  coupling  between  the  deformation  fields  of  the  constituent 
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continua  which  make  up  the  composite.  This  term  is  frequently  neglected 
in  the  study  of  composites^  but  has  been  taken  into  account  in  the  recent 
work  of  McNiven  and  MengiXX  in  their  study  of  two-phase  composites  with 
linear  response. 

Equations  (2.19)  and  (2.20)  together  imply  that 


Ti  =  XiTi  +  XlX2aF(i)E<12^  i=M> 


(2.21) 


where 


T.  : 

i 


Vi}a-  (i)E(i)  <F(i)), 


1,2, 


(2.22) 


are  the  stresses  which  would  arise  at  the  point  X  in  a  single  phase 
medium  composed  solely  of  the  ith  continuum. 


3 ,  Propagation  of  Acceleration  Waves 

In  one  dimension,  the  motion  of  a  nonmaterial  surface  of  discon¬ 
tinuity  with  respect  to  the  reference  coordinates  is  given  by 

Z-Z(t),  (3.1) 

where  Z(t)  denotes  the  position  of  the  surface  in  the  reference  configu¬ 
ration  at  time  t.  The  intrinsic  velocity  U  of  the  surface  of  discontin¬ 
uity  is  given  by 


U(t)  =  >  0  ,  (3.2) 

and  this  quantity  is  a  measure  of  the  speed  of  propagation  of  the  dis¬ 
continuity  surface  with  respect  to  the  reference  coordinates  of  material 
points. 


We  use  the  standard  notation  to  denote  the  jump  in  the  magnitude 
cf  a  quantity  across  the  propagating  surface  of  discontinuity;  thus,  if 
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cp  (X,  t)  is  a  quantity  which  suffers  a  jump  discontinuity  at  the  surface 

<v 

Z  =  Z(t)  but  is  a  continuous  function  of  (X,  t)  jointly  elsewhere,  we 
define  the  jump  in  cp  at  time  t  across  the  propagating  surface  of  discon¬ 
tinuity  to  be 


[cp]  =  [cp]  (t)  =  cp  -cp+  (3.3) 

where 

cp  +  =  Ijim  .  cp (X,  t)  .  (3.4) 

X-*Z  (t) 


Since  U(t)  >  0,  cp*”  and  cp+,  respectively,  denote  the  limiting  values  of  cp 
immediately  behind  and  just  in  front  of  the  propagating  surface.  Of 

12 

course  tep]  must  also  obey  the  kinematical  condition  of  compatibility  , 


(3.5) 


Furthermore,  wc  note  the  formula 

[cpi|t]  =cp+(*l  +  t+£cp]  +  [cpM)]_  .  (3.6) 

A  propagating  nonmaterial  surface  of  discontinuity  is  called  an 
acceleration  wave  if  y(X,t)  and  (X,t)  or,  equivalently,  y^(X,t) 

and  (X,t)  have  the  properties  that  while  v(.,.),  w^  (.,.),  y(.,.), 
w(1)  (.,.),  F ( . , . )  and  F(1)  (.,,)  or,  equivalently,  y(1) (.,.),  y<2]  (.,.), 
y  ^  (.,.),  y(2)  {.,.) ,  F  ^  ( . ,  . )  and  F*2*  (.,.)  are  continuous  everywhere, 
the  second  and  higher  order  partial  derivatives  of  the  fields  y(.,,) 
and  w(1)(.,.)  or,  equivalently,  y^(.,,)  and  y^(.,.)  suffer  jump  dis¬ 
continuities  across  the  propagating  surface  Z=Z(t),  but  arc  continuous 
functions  of  X  and  t  everywhere  else.  Thus,  at  an  acceleration  wave,  we 
have 


ty]  =  tw(1)J  =  [F]  =  [FU)]  -  =  [yU)]  =  [P(1)]  =  [F(2)]  =  [F(2)  J  =  o  . 


(3,7) 
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The  integral  forms  (2.10),  (2.11)  of  the  equations  of  motion  imply 

that  for  all  X  /  2(t)  we  have 

3  k  =  p  y  >  (3*8) 

X  o 

..(1) 

3/  +  5  =  rpow  , 

while  aero? J  the  surface  of  discontinuity  we  nave 


[K]_  +  poU[fr  =  0  , 

(3.9) 

(1) 

L&]  +  rn  U  [w  1=0. 

„  ~  ro  ~  ~ 

It  follows  from  (2.15),  (3.7)  and  the  assumed  continuity  of  the  response 

functions  that  Eqs. (3.9)  arc  satisfied  identically  at  an  acceleration 
wave . 

When  the  jump*  across  the  wavefront  in  Eqs.  (3.8)  are  evaluated  we 


have 

IVL=po^i« 

WjfW  -  rpo[w(1>]^  . 


(3.10) 


It  follows  from  (2.16)  that 


where 


V,-"r2axp*W(1,+V<1>- 


(3.11) 


=  X^*  X2E2  +  ?t1^2(ail+2al2  +  a22>  » 

a2  =  3F(l)K(F,F(1),w(1)  =3/(F,F(1),w(1))=po3F3F(l)2:(B\.F(1),w(1>) 
=  X1EI-rX2E2+X1l?{a11+  (l-r)a12-ra22j  , 
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a3  =  3w(1)^(FjF(1)>w(1))  =- V(f,f(1),w(1))  =  poaFaw(1)r(P,F(1).w(1)) 

=  XxX2  (a13  +  a23)  , 

P2=  9p(l)£(F,F<L),»;(1/)  =PoaF(l)i:(F,F(1)J-w(1)) 


(3.12) 


"  Vl  +  V%  +  X1X2 (ail  '  2rai2  +  r2a22)  ' 


93  =  3wav*<F,FU)jW(1))  --ap(i)5(p,p(1),w(1))  -poap(i)j-w(i)2:(p,Fa),wr-1)) 

-  >./2^12-ra23) 


with 


and 


QT. 

F .  =  — r-rr  ,  no  sum  on  i  , 

1  aF(l) 


(3.13) 


a . 


2— (121 

_a  s  -d) -(2).,  (1), 


53 


(F 


w  )  ,  i,  j  =  1,  2  , 


(3.14) 


which  appear  in  the  last  line  of  each  equation  in  (3.12)  are  in  terms 
of  the  aforementioned  fully  equivalent  alternate  representation.  Since 
the  coefficients  in  (3.12)  are  continuous  functions  for  all  X  and  t,  it 
follows  from  (3.11)  with  the  aid  of  (3.7)  that 

[A  Kl  =<y+  [?I  Fl  +  aV~F(1)]  . 

X  ~  J.i’-’  X  ~  iitv  A  *■** 

(3.3.5) 

I VI = Vvi +  p2l v <1}  l  • 

If  we  put  <p  =  y,  F,  w(1)  and  F^  successively  in  (3.5),  we  find  that 

a  =  r-  ■  Vl =  1,2 IV1 » 

b=  [w(1)^  =-»tF(1)]=U2[Y(l)]  .  (3.16) 


We  call  a  (t)  =  ly]  the  mean  amplitude  of  the  acceleration  x/ave  arid  say  that 
the  wave  is  compressive  if  a  >  0,  expansive  if  a  <  0,  Furthermore,  from 
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(2,3),  (2.8)  and  (3.16)  we  note  that 


!'y(1)]  =  a  +  b,  [y(2)]  =  a  -  rb 


(3.17) 


are  the  amplitudes  of  the  waves  experienced  by  the  two  identifiable 
continua  which  make  up  the  composite. 

The  substitution  of  (3,15)  and  (3.16)  into  (3.10)  results  in  the 
pair  of  coupled  equations 


(p^U2  -  a^)a  -  =  0,  -  r  ’'ay1  +  (p^U2  -  r  1p2>h  =  0, 


(3.18) 


which  admit,  a  nontrivial  solution  in  which  a  f  0,  b^O  provided  U  is  a 
root  of  the  equation 


U4  -  (C2  +  c2)u2  +  (C2C2  -  p)  =  0  , 


(3.19) 


where 


<-r-r  '2=rKV<W>+> 

Ko  o  o  o 


1  +2  1  -i+2  1  ^  +  2 

P  -  ~2~  to2)  =  (d/  )  =  ~  (dF(l)K  )  , 


ind  we  assume  that  >  C xne  roots  of  (3.19)  are 


(3.20) 


U2  =  ~  {  (C2  +  c2)  ±  Jv\  -C2)2  +  43} 


(3.21) 


and,  since  it  is  clear  from  (3,20)2  that  3  >  0,  both  of  the  roots  (3.21) 
will  be  real.  Furthermore, if  we  assume  that 


2  2 

ClS  >  P  ’ 


(3.22) 


then  (3.19)  implies  the  existence  of  two  types  of  acceleration  waves 


whose  speeds  of  propagation  U  ,  U  are  given  by 

r  o 
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UF  C(C1+C2)  +J(Cl-Cl)2  +  4P~3,  (3.23) 

and 

Us  =  1  t  (C1+C2)  -  ^-^)2  +  4P}  ,  (3.24) 

respectively.  Since  C  >  C2,it  follows  at  once  from  (3.22)  -  (3.24)  that 
at  any  point  X  and  time  t. 


U F  ^  C1  and  US^C2,  (3.25) 

with  the  equalities  holding  when  |3=0.  The  term  is  the  speed  one  would 
calculate  from  the  initial  slope  of  the  K-F  curve,  <*  ,  and  thus  (3.23) 
suggests  that  the  " fast"  wave  is  predominantly  associated  with  the  mean 
elasticity  of  the  composite.  Tae  fact  that  U^,  may  exceed  is  a  direct 
consequence  of  the  nonlinear  coupling  effects  which  arise  when  p  does 
not  vanish,  on  the  other  hand,  the  "slow"  wave  always  propagates  into 
a  deforming  composite  behind  the  'fast"  wave  and,  since  (3  is  the  slope 
of  the  curve,  the  "slow"  wave  is  associated  with  the  relative 

motion  of  the  constituents. 

it  has  been  pointed  out  sjy  m.uu&iaeo  and  Walsh  ,  J.ii  a  Somewhat 
different  context,  that  the  inequality  (3.22)  is  capable  of  misrepresenta¬ 
tion.  In  order  that  the  physical  significance  of  the  inequality  may  be 
more  fully  appreciated  we  note  from  Eqs.  (2,12),  (2,13),  (2.16),  (2.19)  - 


(2.22)  and  (3.20)  that 


2  r  2  1  2 

Ci  =  T+7  V1  +  1+7  V2  + 


(l  +  r)p^2) 


{an+2ai2  +  a22], 


(3.26) 


c2  =  -i- 
2  1  +  r 


2  r 

V1  +TTT 


v„  + 


(l  +  r)p(1) 
o 
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},  (3.27) 


and 


15. 


where 


Vi  -d)  V2  ”  “  IT37  sp(2)T2 ’ 


(3.28) 


are  the  intrinsic  velocities  which  acceleration  waves  would  have  in 
bodies  composed  solely  of  either  of  the  interpenetrating  continua  which 
make  up  the  composite. 

It  follows  from  (3.18)  and  (3.20)  that 

b  =  Ha ,  (3.29) 


where 


2  2 

n  “  9  ">  ! 

(U  -C^)r 


(3.30) 


and  it  is  to  be  noted  that  in  the  "  fast"  wave,  (for  which  U  =  U^) , 

H  =  and  sgn  H^sgn  a*  while  in  the  case  of  the  "slow"  wave,  (for 
which  U  =  US),  H  =  Hg  and  sgn  Hs  ^  ~  s9n  &2  •  Equations  (3.17)  and  (3.29) 
together  imply  that 


[y  (1)]^=a(l  +  H)  , 
ty(2)]  =  a (1  -  rH)  , 


(3.31) 


while 


a*>  0  ~  |rya,ip|  >  |[y,2>).pl  , 


[yu,)J  <  |[y‘2)g 


‘“jpl  < 

hy<2llPl , 

u:y  > 

1  ■ 

(3.32) 


(3.33) 


This  is  precisely  the  situation  which  will  always  arise  in  theories  of 
the  type  developed  by  Bedford  and  Stern*5’-^ 

Finally,  let  us  turn  our  attention  to  the  case  of  a  fiber  reinforced 
composite  in  which  the  fiber  is  not  continuous  (i.e.,  chopped  fiber). 
Suppose  that  continuum  2  represents  the  chopped  fiber  continuum.  For  a 
composite  of  this  nature  t2 *=  0  and  consequently  Eq. (2.19)2  implies  that 
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E  =  E(F(1),w(1))  ,  (4.5) 

which,  with  (2.12)  and  (2.19'jj  yields 

K=Po3p(l)2<F(1),w(1))  .  (4.6) 


Under  this  restrictive  circumstance  Eq.  (4,6),  with  (2.12),  (2.13)  and 

(3.20),  leads  to  the  relations 


Cl  =  rC2  =  3F(1^(p(1)’W<1))  ’ 


e  =-2~  (^(i)K)2  =  c^ 


(4.7) 


P  r 
ro 

so  that  Eq. (3.19)  has  only  one  root 

02  =  cJ  +  C2=™cJ=(l  +  r)Cj.  (4.8) 

1  2  r  1  2 

Thus,  in  this  restrictive  case  only  one  wave  propagates  and,  since  H =  1/r 

.•  (2) 

for  this  wave,  from  (3,31)  [y  1=0,  so  that  the  acceleration  of  the 
chopped  fiber  is  continuous  at  the  wavefront. 


5 ,  Variation  of  the  Amplitudes  of  Acceleration  Waves 

In  this  section  we  derive  the  equations  which  govern  the  evolu¬ 
tionary  behavior  of  the  amplitudes  of  acceleration  waves  as  they  propa¬ 
gate  in  two-constituent  composite  materials  modeled  as  interpenetrating 
solid  continua.  We  shall  suppose  that  at  each  instant  both  "fast"  and 
"slow"  acceleration  waves  may  exist  in  the  body.  For  the  moment  there 
is  no  need  for  us  to  distinguish  between  the  two  types  of  acceleration 
waves  nor  do  we  need  to  prescribe  in  detail  the  conditions  which  prevail 
ahead  of  the  waves. 
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The  jumps  in  the  material  time  derivatives  of  Eqs.(3,8)  across  the 
acceleration  wave  yield 

m -rpor«(1)i  .  (5.i: 

On  setting  cp  =  y,  F,  and  successively,  in  (3.5)  we  obtain  the 

relations 

2  ||  -  ~  (Wb+p(1)]-u2pxP(i)i  ,  (5.2: 


the  substitution  of  which  in  (5.1)  yields 


2  ~  =  ~  (  2nV)a  +  —  [Svk]  -  U2  [ftvF]  , 
u  t,  uc  u  ~  X  ~  ~  7 

2§|.~  tfcu)b+7i-  <JVl*H>-o2IV(1,l- 


which  represent  a  set  of  coupled  differential  equations  for  a(t)  and  b(t) 
which  hold  for  each  admissible  propagating  acceleration  wave. 


In  order  to  further  simplify  Eqs.  (5.3)  we  need  to  evaluate  the 
quantities  [i^K]  ,  arid  [-3]  .  Differentiating  Eqs.  (3.11)  with  respect 

to  time  and  evaluating  the  jumps  across  the  wavefront,  with  the  aid  of 
(3.6),  (3.12)  and  (3.16),  we  obtain 

2+  Ta  +?  <s2-a3,b 
1  /■  +  2  +  +  2 

-  —  {c^a  +  2Qfi2at3  +  a22b  3  <5*4> 


?ASa2^  +  ^V  I  flfa  +U  (V2~P3)b 

1  c  +  2  „  +  +  2. 

~  JJJ  W2a  +  2Qf22ab  +  ^22b  3  > 


where 


and 


F(1)\ 

F(1) 

-  "a,)l 

-  ~) 

u  n 

F(1)> 

1  +  a23^ 

VX) 

-  ^1 

~  U  ) 

\ 

u  n 

F{1)S 
U  J 

!  +  M 

-=^)} 

(5.6) 


(5.7) 


with 


<'ii‘ari<r'F<1)'”<1,)  °po''pS:,E''F<1>’"a,> 

■Vl*1!1!*1*  "ill  *  3a112  +  3a122  +  a222'  - 


,(1)  (1) 


(1)  (1). 


.12-3paF(l)K(F,FV  ')-P0Vp(DS(F»P  ’W  5 


-  h\  ~  lX2*2  +  X1X2  (aill  +  (2  "  r)aU2  +  (1  -  2r)a!22  -  ra222)  > 
®13  =  W(F'F<1),tf(1))  ‘pA(^(P'F<1)'W<1)) 


a22 


-  Xj_X2  (a113  +  2al23  +  a223)  , 

^(L)k(F,F(1),w(1))  -  P0^(l)3FZ(F,F(1),w(1)) 


(5.8) 


=  Xl^l  +  r  X2®2  +  X1X2  (alU  +  (1  "  2r)aH2  ”  r  (2  “  r)a122  +  r  a2225  ’ 
a„  =  3c,(l)\,(l)K(F;F(1),w(1))  =pn^(i)Bw(l)BpE(F,F<1),w(1’) 


XlX2(ail3  +  (1  "  r)al23  “  ra223)  ’ 


^22=^(1)^(F’f(1)’w(1))  'P0V1>£(F’P<  >W  > 


^  ^  2  3 

XlFl~r  X2F2  +  XlX2<aill~3ran2  +  3r  al22  **  r  a2.225  > 


222' 


P23=\(1)Bw(1)^(P»f(1)>w(]))  =  Po3F(1>Bw<1>X:(F,P'i,'wVA’> 


.(l)  .  (1), 


while 


XlX2(all3  "  2rai23  +  r  a223  ^  ’ 


32T. 

~  i 

Ei=3i^’ 


dF 


(no  sum  on  i) 


(5.9) 


aV12) 

1^k  G>F(l)SF^)9P<k) 


(F(1)  jF121  ;w(1)  )  ,  i,j,k=l,2 


3a.  . 

a  =  _ LL 

iJ3  .  (1)  7 
3w 


(5.10) 


(5,11) 


which  appear  in  the  last  line  of  each  equation  in  (5.8)  are  :n  terms 
of  the  fully  equivalent  alternate  representation.  It  follows  from 
Eqs.(2.15)3,  (2.16)2  and  (2.18)  together  with  the  definitions  in  (3.12) 

that  we  have 


7=-(a3-wv  j3fS)F-(P3-w'  3p(l.)9)F  +  <SW(1)' 

A  .  UK  .  (1)  -  .  .  UK  '>,••(1) 

w  SwU)9)w  +  <9  +  w  3^U)9)w  , 


(5.12) 


the  jump  in  which,  with  the  aid  of  (3.16),  yields 

rrv-l  +  .  UK  A  +  a  ,rt+  ,  (1)  A+  1)  A  +  ,  (1)  <"•  + 

0]  -  (0f3  - w  3pS  )  -+  (p3~*  (i) 9  )  ~  +  (g  +w  3^(1)9  >b 


(5.13) 


When  the  expressions  (5,4),  (5,5)  and  (5.13)  are  substituted  in 
Eqs.  (5.3)  we  arrive  at  the  coupled  differential  equations 

2H-{rU&0)  ♦rtHro  <K>»-  -7  tv 


-  V*ilV  Icru2,iVi+  r 


(5.14) 


2  f|  =  {§t  +  TTu  (V2'*(i)^U)9+)  +  (g+  +  w(1)3^U)g+)}b 

Mo 

1  r -+  +  •  (1)^  -5l+  -i  1  rr,+  ,2  +  +  2^ 

+  tTv  ^2-a3-W  aFg  > - 3  ^22b  +Vb  +  V  3 

po  rp  U 

+  5“iVl+(c2-u2)lV(1)L-  (5-15) 


■*'“  fyi-sa 
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Of  course,  b  and  a  are  related  through  Eq. (3.29)  and  the  intrinsic 
velocity  U  is  a  root  of  Eq.  (3.19).  l*hese  facts  enable  us  to  combine 
(3.29)  with  Egs.(5.14)  and  (5.15)  in  order  to  obtain  a  single  first  order 
differential  equation  which  governs  the  amplitude  a(t)  of  each  admissible 
acceleration  wave.  After  some  tedious  algebra,  we  find  that  the  amplitude 
of  an  admissible  type  of  acceleration  wave  in  a  two-component  composite 
medium  satisfies  the  equation 

+  pa-  £a2  =  0,  (5.16) 

at 

where 


2(1  +rH?')u  (t)  =  -{  (1  +  rH2)  ~  (toU)+2rH2^  (  foiH)  +  (§*  +  2H|*  + 

-Hw(1)  Opg+  +  HSF{i)g4')/poU-rH2(g+  +  w(1)B^(1)g+)}  (5.17) 


~  3  > 


(5.18) 


2poU 


and 


2(1  +  rh2)  E  =  {a;L1  +  3na12  +  3H2a22  +  H3P22]  =  (1  +  8)1^+  (1  -  rH) 3 

+  i1X2(  (H  H)  3a1L1  +  3  (1  +  H)  2  (1  -  rH)  a112  +  3  (1  +  H)  (1  -  rH)  2a122 

+  (l-rH)3a222]  (5.19) 


is  the  effective  second  order  elastic  modulus  of  the  composite  for  the 
particular  wave  under  consideration,  and  where  the  expression  after  the 
second  equals  sign  is  in  terms  of  the  fully  equivalent  alternate  repre¬ 
sentation. 

Equation  (5,16)  is  a  differential  equation  of  Bernoulli  type.  As 

one  might  expect,  it  is  similar  to  the  equation  recently  derived  by 
16 

Nunziato  and  Walsh  in  their  study  of  the  propagation  of  acceleration 


2 


n 

. 


waves  in  granular  media.  However,  a  cursory  examination  of  the  coeffi¬ 
cients  p,  (t)  and  £(t)  shows  that  the  similarity  is  somewhat  superficial. 

Tlie  coefficients  p.  (t)  and  £(t)  are  determined  by  the  particular  type  of 
wave  under  study,  tire  mechanical  properties  of  the  composite  and  by  the 
conditions  prevailing  ahead  of  the  wave.  We  snail  study  the  properties 
of  the  solutions  of  Eq.  (5.16)  in  a  number  of  particular  situations  in 
the  following  section. 

6 .  The  Behavior  of  Some  Particular  Acceleration  Waves 

in  this  section  we  study  the  evolutionary  behavior  of  the  ampli¬ 
tudes  of  some  particular  acceleration  waves.  In  general,  at  a  given 
instant  of  time,  two  acceleration  waves  will  propagate  in  the  body. 

Suppose  that  the  "fast"  wave  is  located  at  the  poin  X  =  Z  (t)  while  the 

£ 

"slow"  wave  is  at  X  =  Z  (t)  where  Z  (t)  >  Z  (t) ,  In  order  to  simplify 

O  £  s? 

matters  we  shall  assume  that  the  material  ahead  of  the  "fast"  wave  is  in 

a  steady  state  of  equilibrium.  Even  though  the  deformation  field  behind 

a  fast  acceleration  wave  may  be  such  that  (t)  >  U^(t),  the  "slow"  wave 

20  16 

can  never  pass  through  the  "fast"  wave  (cf.  Nunziato  and  Walsh  ). 

Thus,  at  all  points  X  >  Z^(t)  the  fields  y(X)  and  w^(X),  or  alternatively, 
y^  (X)  and  y ^  (X)  do  not  depend  on  t,  i.c., 

y  =  y(x),  W{1)  -  w(1)  (X),  x  >  zp(t)  , 
y(15  -y(1)  (X),  y2  =  y  (2>  (X),  X  >  Z^t)  .  (6.1) 

Since  the  fields  (6.1)  must  satisfy  the  equations  of  equilibrium,  at  all 

points  X  >  Z  (t)  we  have 
F 

"2V  +  W'll,+*I3F<1>-4- 


(6.2) 
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From  (5.16)  we  see  that  the  equation  satisfied  by  the  amplitude  a^tt) 
of  a  " fast"  wave  propagating  into  a  region  which  is  in  a  steady  state  of 
equilibrium  is 

~  aF+^o<t>aF"Co(t)a^=0  ,  (6.3) 


where,  from  (5.37),  with  (5.6),  (5.7)  and  the  fact  that  from  (3.20)  and 
(3.23)  now  U (X)  is  independent  of  t,  wo  have 


2<i  +  *4)Vc  (t)  =-(1  +  rlV>  ¥f  "  "'VWf  +  *1 V 


+  cp2sxp(1)  +cp3f(1)  +rV(0) 


with 


(6.4) 


+2vS)+H?ffS))/pouP' 


^2  (q12  4  2HFa22  +HF^22  )/,poUF’ 


.(0) 


(0)  2Q(°) 


(C+2H„*  '+«  P  V)/pA, 


__  ,  2  (0)w  (0) 
HF  =  (PqUF  ‘  “l  )/a2  > 


(6.5) 


and  C  (t)  is  still  given  by  (5.18).  The  superscript  0  occurrina  in  g^ 
o 

and  on  the  right-hand  sides  of  the  expressions  in  (6.5)  denote  that  these 

quantities  arc  evaluated  for  the  steady-state  deformation  fields  described 

by  (6.1)  and,  consequently,  these  quantities  are  functions  of  X  only  and 

do  not  depend  on  t.  Equation  (6.3)  has  the  same  form  as  that  which 

governs  the  evolutionary  behavior  of  the  amplitude  of  one-dimensional 

acceleration  waves  in  a  single  phase  clastic  material  which  is  in  a  state 

_  12 

of  nonhomogrneous  deformation  ahead  of  the  wave  (see,  e.g.,  Chen  , 

.  21 

Coleman,  Greenberg  and  Gurtxn  ) . 

If  the  material  ahead  of  the  "fast"  wave  is  at  rest  in  its  natural 


stress-free  state  so  that 
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y  =  x,  v(1)  =  o,  x  >  zp(t)  , 

or 

Yi(x)  =y2(x)  =  x »  x  >  zF(t)  * 
and  it  follows  from  (6.4)  that 


2  (0)  2 

V-Q  <t)  =  rHpg  /2  (1  +  rHp)  =  u)o  ,  say, 

for  such  a  wove  so  that  a  (t)  obeys  the  differential  equation 

F 


2 

— — -  =  -  cu  a  +  C  a  , 
at  o  F  F  > 


(6, 


(6.8) 


(6.9) 


The  properties  of  the  solution  (•'  .10)  are  well  documented  (see, 
e.g..  Refs. 12  and  23)  and  it  is  not  our  intention  to  study  them  in  detail 
here.  The  critical  mean  amj>litude  for  acceleration  waves,  Xq,  plays  a 
fundamental  role  in  determining  whether  the  amplitude  of  an  acceleration 
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wave  will  grow  or  decay  as  the  wave  traverses  the  material.  In  particular, 
the  sign  of  XQ  plays  a  critical  role  in  determining  whether  the  mean 
amplitude  will  grow  or  decay.  Since  g^  >  0,  by  assumption,  it  follows 
from  (6.8)  that  >  0  so  that  sgn  (X  )  =  -  sgn  (E)  and  it  is  clear  from 
(5,8)  and  (5.19)  that  E  may  differ  in  sign  from  either  or  E^.  Notice 
that  (6.10)  implies  that 

(i)  If  J a^ ( 0 )  }  <Xq,  then  a  (t)  —  0  monotonically  as  t-><=. 

(ii)  If  sgn  (a^  (0)  )  =  sgn  and  j  ( C )  1  >  |  X  [,  then  (t)  -*  oo 
monotonically  within  a  finite  time 

to  =-  ~  lull-  (Xo/aF(0))}  ,  (6.12) 

o 

and  this  is  usually  taken  to  indicate  shock  formation. 

Next,  let  us  consider  an  acceleration  wave  propagating  into  a  com¬ 
posite  material  in  which  3  is  independent  of  w*1^  so  that  g(0)  =0. 

Equation  (6.9)  now  reduces  to 


at 


2p 

F 


3  aF 


=  0 


(6.13) 


T.T  V>  *C  —I  flvn  ■14-*-’  4-V.  r  -  1  %  1  4-  *»  / 


aF(t) 


aF(0) 


(6.14) 


1  + 


a  (0)E 
F _ 

3 

2o  DT 
’  o  F 


If  sgn  (a  /q)E)  <  0,  then  the  solution  (6,14.)  becomes  unbounded  after  a 


time 


-2p  IT 

t=  -J°-£ 

aF(0)E 


(6.15) 


V,  1 
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and,  of  course,  this  is  precisely  what  happens  in  a  perfectly  elastic 
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single  phase  continuum  (see,  e.g.,  Green  ). 

If  the  material  ahead  of  the  fast  wave  is  at  rest  such  that  the 
center  of  mass  y  is  in  a  state  of  homogeneous  strain,  we  have 


y=XX,  X  =  constant,  X>ZF(t)  .  (6.16) 

Since  F  =  X,  a  constant,  in  this  case  it  follows  from  Eqs,  (6.2)  that  at 

all  points  X>Z  (t)  the  relative  displacement  w^  (X)  must  satisfy  the 
F 

nonlinear  first-order  differential  equation 


dw<D 

(a332_33a2)  dx 


(6.17) 


and  it  is  to  be  noted  that  both  and  the  coefficients  occurring  in  (6.17) 
are  (for  fixed  l‘=X),  functions  of  and  w^.  It  is  to  be 

expected  that,  even  if  g+(0)=0,  the  solution  of  (6.17)  for  w*]^(X)  will 
load  to  a  nonvanishing  expression  for  pQ(t)  when  substituted  into  (6.4). 
Thus,  the  coefficient  ^(t)  will  be  a  consequence  of  the  inhomogeneities 
in  the  deformation  fields  of  the  two  continua  which  make  up  the  composite, 
as  well  as  the  coefficient  g+(0).  The  solution  of  Eq.  (6.3)  in  thi.s 


The  properties  of  the  solution  have  been  discussed  in  detail  by  a  number 
of  authors  (Bailey  end  Chen  ,  Nunziato  and  Walsh  )  and  we  refer  the 
reader  to  these  works  for  details. 

The  behavior  of  "slow"  waves  is  always  more  complicated  than  that 
of  "fast"  wavec  since  "slow1'  waves  propagate  into  regions  which  are  not 
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in  equilibrium.  Thus,  in  the  case  of  the  "slow*'  wave  the  coefficients 
^ (t)  and  £(t)  will  contain  derivatives  with  respect  to  both  X  and  t. 

The  amplitude  of  the  "slow"  wave  is  given  by  an  expression  similar 
to  (6.18). 

In  Section  4  we  examined  the  propagation  of  waves  in  a  number  of 
special  situations.  Let  us  now  consider  the  evolutionary  behavior  of 
the  amplitudes  of  these  waves. 

We  saw  in  Section  4  that  in  the  restrictive  case  when  the  internal 
energy  is  given  by  (4.1)  two  waves  may  propagate  with  intrinsic  veloci¬ 
ties  U  =  V  and  U =  V  ,  respectively.  The  equations  satisfied  by  the 

FI  S  2. 

amplitudes  tf  these  waves  may  be  deduced  from  the  results  of  Section  5, 
but  the  properties  of  the  waves  in  the  special  highly  simplified  cases 
considered  here  become  more  transparent  when  it  is  noted  that,  with  the 
aid  of  (2.3),  (2.8),  (2.9),  (2.12),  (2.13),  (2.21)  and  (2.22),  Eqs.  (3.8) 
may  be  transformed  to 


,  „  .  1  „  -(1)..(1) 
3xTi  +  l-uTTT^Po  y  ’ 


-s  m  1  „  _(2)..(2) 

SXT2  x.  fl  +  r)  3  Po  Y 
2 


(6.19) 


It  follows  from  (2.16)2  and  (4.1)  that 

3F--P  X.Ud  (1)E(12)  (w(1))  ,  (6.20) 

o  1  /  w  7 

A  .  (1) 

and  we  further  assume  that  3=  gw  ,  where  g  is  a  negative  constant.  If 
we  write  ly^)  =  a^  and  =a^,  then  an  elementary  calculation 

leads  to  the  growth  equations 


d5(l)  .  ..  (i)_ (i)  „(i'_  (i)2 

“dt“  +  u  a  -C  a  ■ 


0, 


i-1,2  , 


(6.21) 
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where 


2Pq  x 

£(i)=VPoi)vi>  V^V  i==1>2* 


(6.22) 


It  follows  from  (6.19)  ,  that  at  the  "fast"  wave  we  have 


r  3  (2)  3 _ *1* _  (1) 

[3  y  /3t  1,  22  2  a 

F  \JpU)  a  +  T)  iv:-vz_) 

2ro  2.  1 


(6.23) 


so  that  the  third  and  higher  order  derivatives  of  y  (X, t)  suffer  jump 
discontinuities  across  the  "fast"  wave.  Likewise  at  a  "slow"  wave 
y(±>  (X,t)  and  its  first  and  second  partial  derivatives  will  be  continuous 
fojT  all  t j  £\ifc 


oV11/^ - -  -»,2> 

h'o  (1+r> 


(6.24) 


On  the  other  hand,  in  materials  in  which  g  =  0,  the  right-hand  sides  of 
the  expressions  (6.23),  (6.24)  vanish  and  the  fourth  order  time  deriva¬ 
tives  of  y*2^  and  y^,  respectively,  suffer  jumps  across  the  wavefronts 
which  are  given  by  the  expressions 


.  4  (2)  4  ..  Vw(1),? _  , 

[3  y  /3t  ]  -  »2\  2  2  1 

Vo  <1  +  r)<Vl-V2> 


(6.23a) 


=  _ Vw*1^  , 

[3  y  /3t  /ns  5  JT  ' 

s  x.p(1)  (1  +  r)  (V2-V2) 
1  o  2  1 


(6.24a) 


Hie  importance  of  the  foregoing  is  that  in  a  composite  of  the  type 


characterized  by  the  restrictive  expression  (4.1)  for  the  internal  energy 
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per  unit  mass,  a  "fast"  wave,  moving  into  a  material  which  is  initially 

at  rest  in  an  equilibrium  configuration,  will  induce  a  motion  in  both 

constituents  of  the  composite.  Notice  from  (6,22)  that  the  evolutionary 

behavior  of  a  particular  wave  is  influenced  only  by  the  properties  of 

one  of  the  constituents  and  the  state  of  this  constituent  ahead  of  the 

wavefront.  Thus  we  note  that,  in  particular,  the  behavior  of  a  "fast" 

wave  propagating  into  a  composite  which  is  at  rest  in  a  homogeneous  state 

before  the  arrival  of  the  wavefront  is  qualitatively  the  same  as  that  of 

a  wave  propagating  into  an  equilibrium  configuration  in  a  single  phase 

thermoelastic  medium.  On  the  other  hand  the  second  constituent  of  the 

composite  will  be  set  in  motion  through  the  coupling  caused  by  3  because 

(2) 

of  the  passage  of  the  "fast"  wave  and  consequently  |i  will  generally 
be  nonzero  even  when  g  =  0. 

The  expressions  (6.23),  (6.24),  (6.23a),  (6.24a)  are  also  interest¬ 
ing  in  that,  when  the  internal  energy  is  given  by  (4.1).  it  is  evident 
that  the  higher  order  discontinuity  induced  in  the  motion  of  one  com¬ 
ponent  of  the  composite  because  of  a  discontinuity  in  the  acceleration 
of  the  second  component  across  the  wavefront  is  the  result  of  coupling 
effects  caused  by  the  relative  body  force  3. 

To  complete  our  study  we  return  to  the  case  of  a  composite  in  which 
one  of  the  components  is  a  chopped  fiber.  The  internal  energy  density  is 
now  given  by  (4.5)  and  only  one  wave,  across  which  /  0,  ly^]  =0, 

may  propagate  in  the  composite.  When  Eqs.(2.3),  (2.8),  (2.12),  (2.13) 
and  (3.8)  arc  combined,  it  follows  that  the  motion  of  the  chopped  fiber 
(i.e.,  component  2  of  the  composite)  is  given  by  the  formula 


=  -3  , 


(6.25) 


I 


5  =_  P0aw<1)^(^<1)>w(1))  +^<1}  > 


(6.26) 


and  once  again  we  assume  that  g  is  a  positive  constant.  It  follows  from 
(6.25)  and  (6.26)  that 


[dV2)/at3]  =-  —  m  =  fa-(l)^-gu}a  , 


(6.27) 


2  hi  2 

■hhere  a  =  [5  y'  J/Bt  ] . 

The  cne-dimensional  behavior  of  a  eentrosyvnmetric  medium,  of  the 
type  under  study  is  characterized  by  an  internal  energy  density  function 
j,(l)  (j,(l)^w^^ )  which  is  an  even  function  of  w^,  i.e., 


2(P(1),W(1))  =.S(P(1),  -W(1))  , 


(6.28) 


~  “ ^1)  fl) 

so  that  the  relative  body  force  3(FV  ,w'  )  is  an  odd  function  of  w'*'. 

It  follows  that  )3(F^1\o)-(i)_1  =  0  so  that  if  the  material  ahead 
of  the  wave  is  in  equilibrium  in  its  natural  stress-free  state  then 


3  <21  3 

£5  y  Vat  ]_=ga/por  , 


(6.29) 


3  (2)  3 

while  if  g  =  0,  By  /St  is  continuous  for  all  X,t  but 


4  (2)  4 

[B  y '  /at  ]  =  a  (l)!?a/p  r  . 

—  —  w  o 


(6.30) 


Finally,  let  us  consider  the  behavior  of  the  amplitude  of  an  accel¬ 
eration  wave  which  is  propagating  into  a  composite  in  which  one  of  the 
components  is  a  chopped  fiber  and  which  is  in  a  steady  natural  stress- 


free  state  before  the  arrival  of  the  wave.  It  follows  from  Section  5 
that  the  amplitude  of  the  wavefront  is  governed  by  Eq. (5.16)  with 


2  (1  +  r)  q>  3  (1  +  1/r)  (Xl\  +  X2aill)’ 

U 

po 


(6.31) 
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so  that  at  anytime  t  the  amplitude  of  the  wave  is  given  by 


a  (t) 


X 


o 


(6.32) 


where  Xq  =  £T/£.  Thus,  in  a  chopped  fiber  composite  the  behavior  of 
the  acceleration  wave  is  qualitatively  the  same  as  that  in  a  single 
phase  heat  conducting  elastic  medium.  Notice  in  particular  that  the 
wave  will  be  undamped  if  g  =  0,  in  which  case  (6.32)  reduces  to  the 
expression 


aft) 


a  (0) 

1  +  Ca.  (0)  t 


(6.33) 


The  influence  of  the  chopped  fiber  on  the  behavior  of  the  amplitude  is 
evident  from  the  manner  in  which  the  parameters  r,  X2  and  g  influence 
the  coefficients  gl  and  Q. 
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